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Summary 

Long-period variations have been observed in the mean values of ec- 
centricity and argument of perigee (published by the Goddard Space Flight 
Center) for Nimbus 2. It is shown that the use of a more recent determina- 
tion of the earth zonal harmonic coefficients, plus additional first-order 
earth oblateness perturbations (for small eccentricity satellites) not in- 
cluded in the mathematical model for orbit determination, explains these 
variations. 
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LONG-PERIOD VARIATIONS IN THE ECCENTRICITY 
AND ARGUMENT OF PERIGEE OF NIMBUS 2 


INTRODUCTION 

In two earlier papers (References 1 and 2), the problem of small divisors 
in the case of low eccentricity satellites was analyzed at some length. Eccen- 
tricity was found to occur as a divisor in analytic expressions (both long and 
short period) for the perturbations, caused by the non-sphericity of the earth, in 
eccentricity, mean anomaly, and argument of perigee. Formulas for these per- 
turbations were presented, and contain terms involving up to J \ and (J { /J 2 ) 3 , 
where J . is any odd zonal harmonic. In addition, the long-period terms were 
found to explain satisfactorily the observed variations in eccentricity and argu- 
ment of perigee for the satellites Alouette 1 and Tiros 8. 

Similar long-period variations in eccentricity and argument of perigee were 
observed in the case of Nimbus 2, another low eccentricity satellite. It will be 
shown here that the long-period expressions in References 1 and 2 are adequate 
representations of these perturbations. 


LUNAR AND SOLAR EFFECTS 

Mean values for the eccentricity and argument of perigee, developed at the 
Goddard Space Flight Center, were first corrected for lunar and solar gravita- 
tional .effects and the effects caused by solar radiation pressure and lunar and 
solar tides; the formulas used to compute these perturbations appear in Refer- 
ence 3 and 4. For the most part, these effects were fairly small, except for 
three near-resonant solar terms in the argument of perigee perturbations. 
However, since these terms all had periods in excess of 10 4 days, their effects 
over the time interval studied would appear to be secular and therefore wouldn't 
affect the analysis of the periodic variations. 

The corrected values of eccentricity and argument of perigee appear in 
Tables 1 and 2 as e" and g" , respectively. 
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LONG-PERIOD ZONAL HARMONIC PERTURBATIONS 


Analytic formulas for the major long-period zonal harmonic effects on the 
eccentricity and argument of perigee for a low eccentricity satellite appear in 
References 1 and 2; they are reproduced here for convenience: 
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where a", e p and i" are semi-major axis, eccentricity, and inclination con- 
stants, respectively. In addition, 
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g" = argument of perigee 



2 


where 


3 J 2 , 3J2 

N = (1-5 cos 2 i") + i — (7 - 114 cos 2 i" + 395 cos 4 i") 

4a" 3 v / a" 64a" s v / ? r 


15J- 

(3 - 36 cos 2 i" + 49 cos 4 i") + . . . 

32a" s v^ r 


3J, sin i" 

M = - (1 

Sa'^VaT 


, 15J S sin i" 

5 cos 2 i") + (1 - 14 cos 2 i" + 21 cos 4 i") 

32a" 6 v / a 57 * 


105J - sin i" 

+ (5 - 135 cos 2 i" + 495 cos 4 i" - 429 cos 6 i") 

1024a" 8 '/a 7r 


315J g sini" 

+ (7 - 308 cos 2 i" + 2002 cos 4 i" -4004 cos 6 i" + 2431 cos 8 i") 

4096a" 1 % r 


3465J., sini" 

+ — (21 - 1365 cos 2 i" + 13650 cos 4 i" - 46410 cos 6 i" 

131072a" 12 v / a"~ 


+ 62985 cos 8 i " - 29393 cos 10 i") + 


However, the terms -Q cos 6 (for Ae ) and (Q /e x ) sin 6 (for Ag ) have already been 
used in the mathematical model for orbit determination, at least for the zonal 
harmonics J 2 , J 3 , J 4 , and J g . The following values were used for J 3 and J s : 


J 3 = -2.285 x 10" 6 , J 5 = - 0.232 x 10“ 6 . 


It is suggested that Kozai's 1964 determination of the harmonic coefficients (see 
Reference 5) form a better set - they are 
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1.082645 x 10“ 3 

J 3 = -2.546 x 

10“ 6 

-1.649 x 10‘ 6 

J 5 = -0.210 x 

10“ 6 

-0.333 x 10 -6 

J 9 = -0.053 x 

io - 6 

0.302 x 10 -6 . 




It is further suggested that the observed variations in e and g would reflect the 
differences between the more accurate J 3 , J 5 values and the values which were 
actually used in the -Q cos 6 and (Q /e 1 ) sin 0 terms. 

Therefore, setting 

AJ 3 = -2.546 X 1CT 6 - (-2.285 x 10 -6 ) = -0.261 x 10' 6 
AJ = -0.210 x 10' 6 - (-0.232 x 10 -6 ) = 0.22 x 10 -7 
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the expressions to be compared with the observed variations are 
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Using 


a" = 1.1782076 earth radii 


e j = .00560 
i" = 100S306 

and Kozai's 1964 determination of the zonal harmonic coefficients, Equations (4) 
become 

Ae = -(1.092397 lx 10' 4 ) cos 6 - (0.63515970x 10“ 4 )cos 26 + (0.85403736x 10 ~ 5 )cos 38 
Ag = 1?3454747 sin 6 + 1P2997132 sin20 +0^23301303 sin 36 . (5) 


OBSERVED VARIATIONS 

The values of e" and g" were then fit, by means of least squares, to the 
following models: 



A cos n 8 
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where e Q , g 0 , A n , and B n are constants, g". is the secular motion of g" , t - t 0 
is the elapsed time in days, and 


6 = 51?34033 - (2?3536297/day) (t - t Q ). 
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The constant 51°.34033 is merely the first value of g" plus 90°, and the secular 
motion of 8 was computed from the expression for dg" /dt appearing in Refer- 
ence 6. The results are 

e".= 0.55936191 x 10" 2 - (1.200321 x 10“ 4 ) cos 0 - (0.648964 x 10" 4 ) cos 28 
+ (0.44700 x 10' 5 ) cos 38+ (1.5399* 10' 6 ) cos 45+ (1.7144x 10 -6 ) cos 58 
- (0.9668 x 10' 6 ) cos 6 8 + (0.34168 x 10 -5 ) cos 78 +(1.9838 x 10 -6 ) cos 88 
+ (0.8763 x 10 -6 ) cos 98 


g'.= 679P62285 - (2?3542 1 10/day ) (t - t Q ) + 1?3007803 sin 8 
+ 0?8192370 sin 28 - 0?0202026 sin 38 + 0?0 196037 sin 4 8 
+ 0?0099007 sin 58 + 0?0441686 sin 6 8 + 0f0034652 sin 78 
+ 0?0379 263 sin 86 + 0?0279919 sin 98. 


A comparison between Equations (5) and (7) indicates that the major observ- 
able periodic variations in e and g (the 8 and 28 terms) are explainable by the 
theory. There is somewhat of a discrepancy between the amplitudes of the sin 28 
terms. However, it should be noted (see Equations (4)) that the theoretical am- 
plitude is sensitive to the value of e x used; use of a larger value would bring 
closer agreement. 

Figures 1 and 2 show the closeness of the least squares fits. 
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Table 1 


Eccentricity of Nimbus 2 


t - t Q (days) 


e".x 10 



A cosn<9 


10 


0 

.55332889 

7 

.54728658 

14 

.54168522 

21 

.54233060 

28 

.54606277 

35 

.54525536 

42 

.54996346 

49 

.55606657 

55 

.55957010 

62 

.56235129 

69 

.56480941 

76 

.56463475 

83 

.56371704 

90 

.56290166 

97 

.56267565 

104 

.56470174 

111 

.56764970 

118 

.57106100 

125 

.57244757 

132 

.57036198 

139 

.56764492 

146 

.56238443 

153 

.55574878 

160 

.54816344 

167 

.54255936 

174 

.54007409 

181 

.54265645 

188 

.54266910 

195 

.55178840 

202 

.55480257 

209 

.56466682 

216 

.56788385 

223 

.56750003 

230 

.56620046 

237 

.56768515 

244 

.56621362 

251 

.56548537 

258 

.56634020 

264 

.56908258 

271 

.57008731 

278 

.57087450 

285 

.56981932 

292 

.56534614 

299 

.55772803 

306 

.55368415 

,0055936191 

Aj = -.0001200321 


.55309715 

.54689101 

.54263530 

.54216936 

.54230895 

.54553649 

.55172237 

.55742556 

.56257593 

.56728237 

.56868472 

.56837489 

.56742564 

.56565369 

.56413156 

.56492539 

.56686812 

.56810100 

.56866868 

.56814157 

.56465560 

.55874767 

.55306014 

.54685197 

.54262321 

.54216969 

.54231507 

.54557165 

.55176106 

.55746235 

.56344980 

.56768853 

.56869821 

.56827280 

.56721607 

.56534977 

.56408524 

.56520973 

.56687891 

.56810638 

.56867023 

.56813067 

.56462217 

.55870967 

.55302309 

= -.0000648964 A 3 = .0000044700 


A 4 = .0000015399 A s = .0000017144 A 6 = -.0000009668 A ? = .0000034168 

A g = .0000019838 A g = .0000008763 6 = 51^34033 - (2?3536297/day) (t - t 0 ) 

t - t 0 - days since 5/22/66, 0 hrs. U.T. 
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Table 2. Argument of Perigee of Nimbus 2 
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0 (days) 

g"(deg.) 

K ~ [ go + Sc ( t - t o^ ( de S- > 

7 B s in i 

L _■ n 

n - 1 

0 

681.34033 

1.71748 

1.81008 

7 

664.56280 

1.41943 

1.42397 

14 

647.43303 

0.76913 

0.98474 

21 

630.18123 

-0.00319 

0.12833 

28 

612.91141 

-0.79353 

-0.84045 

35 

596.17881 

-1.04665 

-1.32491 

42 

579.25093 

-1.49506 

-1.74587 

49 

562.60975 

-1.65676 

-1.81638 

55 

548.70462 

-1.43662 

-1.60563 

62 

532.58137 

-1.08040 

-1.22130 

69 

516.59107 

-0.59122 

-0.63771 

76 

500.47859 

-0.22422 

-0.25701 

83 

484.09840 

-0.12494 

0.00321 

90 

467.60543 

-0.13843 

0.15816 

97 

450.82827 

-0.43611 

0.03784 

104 

434.13663 

-0.64828 

-0.14006 

111 

417.82345 

-0.48198 

-0.08924 

118 

401.75990 

-0.06605 

0.16244 

125 

385.67054 

0.32406 

0.46239 

132 

369.87252 

1.00552 

1.01109 

139 

354.07590 

1.68838 

1.49166 

146 

337.83500 

1.92696 

1.77110 

153 

321.44596 

2.01739 

1.80866 

160 

305.04411 

2.09502 

1.42115 

167 

287.97163 

1.50202 

0.98117 

174 

270.32661 

0.33647 

0.12129 

181 

252.80550 

-0.70516 

-0.84488 

188 

236.12142 

-0.90976 

-1.32760 

195 

218.76501 

-1.78669 

-1.74801 

202 

202.22859 

-1.84364 

-1.81532 

209 

186.04263 

-1.55012 

-1.56182 

216 

169.99133 

-1.12194 

-1.13988 

223 

153.95519 

-0.67861 

-0.56128 

230 

137.79547 

-0.35885 

-0.21829 

237 

121.45899 

-0.21585 

0.04051 

244 

105.04672 

-0.14865 

0.15577 

251 

88.51341 

-0.20248 

0.00728 

258 

71.94073 

-0.29568 

-0.15223 

264 

57.97534 

-0.13581 

-0.08798 

271 

41.83007 

0.19840 

0.16406 

278 

25.69696 

0.54477 

0.46502 

285 

9.76031 

1.08759 

1.01492 

292 

- 6.29829 

1.50843 

1.49380 

299 

-22.50038 

1.78586 

1.77253 

306 

-38.95179 

1.81392 

1.80722 


g". = -2.°3542110/day 
B . = 0.°0196037 


B 0 = 679?62285 
B 3 = -0.°0202026 
B, = 0?0034652 

6 = 51? 34 033 - (2?3536297/day) (t - t 0 ) 


B g = 0?0379263 


Bj = 1?3007803 B 2 = 0?8192370 

B s = 0?0099007 B g = 0?0441686 

B 9 = 0?0279919 
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